Experimental test of a trace formula for two-dimensional dielectric resonators 
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Resonance spectra of two-dimensional dielectric microwave resonators of circular and square 
shapes have been measured. The deduced length spectra of periodic orbits were analyzed and 
a trace formula for dielectric resonators recently proposed by Bogomolny et al. [Phys. Rev. E 78, 
056202 (2008)] was tested. The observed deviations between the experimental length spectra and 
the predictions of the trace formula are attributed to a large number of missing resonances in the 
measured spectra. We show that by taking into account the systematics of observed and missing re- 
sonances the experimental length spectra are fully understood. In particular, a connection between 
the most long-lived resonances and certain periodic orbits is established experimentally. 



PACS numbers: 05.45.Mt, 42.55.Sa, 03.65.Sq 

I. INTRODUCTION 

Trace formulas relate the density of states of a 
quantum system to the periodic orbits (POs) of the 
corresponding classical system. They were first intro- 
duced by Gutzwiller jlHSj and have since then found 
numerous applications not only to quantum systems but 
also to wave-dynamical systems such as electromagnetic 
m |3] or acoustic [B] resonators. In addition concepts 
from quantum chaos and semiclassics are nowadays 
also applied to open dielectric resonators (also called 
dielectric billiards), which are used, e.g., as microlasers, 
as sensors or in optical circuits [3 H]. Especially the 
occurrence of scarred resonance states [H [10] and the 
connections between the emission properties of micro- 
lasers and POs have gained considerable attention. It 
has been established that the directions of maximal 
emission from a microlaser with, e.g., quadrupole shape 
are determined by the unstable manifolds of certain POs 
[TTl [T^ . Recently, a trace formula for two-dimensional 
(2D) dielectric billiards was proposed in [TJ. This trace 
formula provides a connection between the resonance 
density of a dielectric resonator as a wave-dynamical 
system and the POs of the corresponding classical 
ray-dynamical billiard system and is a continuation of 
ideas developed in [H] and [15]. 

The objective of the present work is an experimental test 
of the trace formula [T^] with 2D dielectric resonators. 
Due to radiation losses, the resonances of a dielectric 
billiard have finite lifetimes. The lifetimes Tj depend 
on the individual resonances, and the corresponding 
resonance widths Tj = 1/tj range over several orders 
of magnitude. Therefore, the spectrum consists of a 
mixture of resonances with very small up to extremely 



large widths, and only the sharp, long-lived resonances 
can be clearly identified. Thus, only a part of the 
resonances is actually observed in an experiment. This 
is an important difference to hardwalled billiards, for 
which complete spectra can be measured using super- 
conducting microwave resonators [S] [THl [TT] . The task of 
the present work is to demonstrate that nevertheless the 
measured, incomplete spectrum can still be interpreted 
by means of the trace formula. Furthermore, long-lived 
resonances, which play an important role in microlasers, 
and their relation to certain POs are investigated. It 
should be noted that flat microlasers are usually ap- 
proximated as 2D systems by introducing the so-called 
effective index of refraction, even though the precision 
of this approximation is not always under control [18] . 
The applicability of the trace formula on such flat 
three-dimensional systems within this 2D approximation 
will be the subject of a future publication. 
Microwave resonators provide a suitable testbed for the 
investigation of dielectric resonators because of their 
macroscopic dimensions and the large spectral range 
accessible experimentally. Here we use flat dielectric 
plates of different shapes and materials as passive 
resonators and put them between two copper plates. 
Then, up to a certain frequency, they are described by 
a 2D scalar Helmholtz equation [16]. Three different 
resonators with regular classical dynamics were inves- 
tigated, a circular and a square one made of Teflon 
each with index of refraction n « 1.4 and a square 
resonator made of alumina (AI2O3) with n w 3. The 
paper is organized as follows. Details of the trace 
formula are briefly summarized in Sec. |IT| and the 
experimental setup is described in Sec. IIIII The results 
for the circular Teflon resonator, the square Teflon 
resonator, and the square alumina resonator are pre- 
sented in Sees. [IV[ |V] and [VI[ respectively. Section |VII| 
concludes with a discussion of the results and a summary. 
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II. TRACE FORMULA FOR DIELECTRIC 
RESONATORS 

As already noted above, trace formulas relate the den- 
sity of states of a wave-dynamical system to a sum over 
all POs of the corresponding classical ray-dynamical sys- 
tem. An open 2D dielectric resonator is a flat cylinder 
whose cross-sectional area has an arbitrary shape made 
of a dielectric material with index of refraction n > 1 
surrounded by air (or another material with lower in- 
dex of refraction). The corresponding classical system 
is the billiard with the same shape. Rays travel freely 
inside the billiard domain S and are partially reflected 
back inside and partially transmitted outside of the bil- 
liard according to the Fresnel formulas [TH] when they hit 
the boundary dS. The wave equations used for such 2D 
dielectric resonators are [201 
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The wave number k is related to the frequency / via 
k = 27r/ / c, where c is the speed of light in vacuum. Be- 
low a certain frequency only transverse magnetic (TM) 
fleld modes are excited in the resonator setup considered 
here (see Sec. IIII|) . Accordingly, the wave function 'i/ 
corresponds to the ^-component of the electric field, E^, 
and both the wave function and its normal derivative 
are continuous along the boundary dS. In such a res- 
onator there are only quasi-bound states or resonances, 
which are characterized by complex frequencies fj , where 
Re(/j) is the resonance frequency and Tj = — 2Im(/j) is 
the resonance width (full width at half maximum) [7]. 
The resonance frequencies and widths are obtained by 
fitting Lorentzians to the measured frequency spectra. 
The spectral density of states is defined as 
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The density of states can generally be decomposed into a 
smooth part (also known as Weyl term [5T|) and a fluc- 
tuating part, p — pvvcyi + Pfluc, where the smooth part is 
related to the area A and perimeter U of the resonator 
and the fluctuating part to the POs of the classical bil- 
liard. For a 2D dielectric resonator the smooth part of 
the density of states is [13. 
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Here, f is related to the boundary conditions. It is de- 
fined as 



f{n) = 1 + 
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where for TM modes 
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In the case of a regular classical billiard the semiclassical 
expression for pfluc reads as |13j 
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Here, 



Bpo oc 



where Apo 



is the area of the 



billiard in configuration space covered by the family of 



POs with length 



Rpo denotes the product of all 



Fresnel refiection coefficients resulting from reflections at 
the boundary of the billiard, and ippo denotes the phase 
changes accumulated at the reflections [i.e., arg(i?po)] 
and the caustic points [131 [21]. The details of £po, Bpo, 
ippo and Rpo for the dielectric circle and square billiards 
are given in Appendixes [X] and [C] For 7~ = — 1 this trace 
formula equals that of a closed 2D quantum billiard with 
Dirichlet boundary conditions except for some additional 
terms involving the index of refraction n which reflect the 
larger optical length or volume of the resonator and the 
contributions of the Fresnel coefficients. The semiclassi- 
cal expression pweyi + pIuc ^ good approximation to 
the density of states p{k) defined by Eq. ^ in the semi- 
classical limit fc — >■ oo. The Fourier transform (FT) of 

Pfluc(fc), 
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yields the length spectrum which has peaks at the lengths 
£po of the periodic orbits. For a test of the trace formula 
we compare it with the FT of the fluctuating part of the 
density of states. 
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The summation in Eq. (jsj is over all resonances with 
fcmin < ^^(fcj) < fcmax- Short-livcd rcsouanccs are sup- 
pressed because of the factor exp (— n/|Im(fcj)|) appearing 
in Eq. (|8| as compared to long-lived resonances. 



III. EXPERIMENTAL SETUP 

A sketch of the experimental setup is shown in Fig. [T] 
a plate made of a dielectric material (Teflon or alumina) 
is placed between two copper plates. Then below the 
frequency 



/2D — 
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(9) 



only TM modes with homogeneous electric field in the 
0-direction (TMq) exist. Here, d is the thickness and n 
is the index of refraction of the dielectric plate. Below 
this frequency, the resonator can be treated exactly as a 
2D system [16 . Microwave power is coupled into and out 
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(a) Side view: the dielectric plate is placed between 
Two dipole antennas entering the res- 



FIG. 1. Schematic picture of the experimental setup (not to 
scale) 

two copper plates, 
onator through small holes in the top plate are placed next 
to the sidewalls of the dielectric plate. The dipole antennas 
are attached to a vectorial network analyzer via rf cables. |(b)| 
Top view: the solid lines denote the contour of the copper 
plates; the dashed lines denote that of the dielectric plates 
(circle with radius R and squares with side length a). The 
crosses indicate the positions of the antennas. 



of the resonator with two dipole antennas. A vectorial 
network analyzer (PNA 5230A by Agilent Technologies) 
is used to measure the complex transmission amplitude 
'S'21 (/), where the modulus squared of S21 equals the ratio 
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between the power Pjn coupled in by antenna 1 and the 
power Pout coupled out via antenna 2 for a given fre- 
quency /. The dipole antennas are put next to the 
sidewalls of the dielectric plates [see Fig. l(a)| , so that 
they can couple to the evanescent fields of the resonance 
modes. Their positions are indicated in Fig. 1(b) The 
measured transmission amplitude 15211 at a resonance 
frequency rises with the electric field strengths at the 
positions of the two antennas [22]. In the case of the 
circle billiard, oc cos (rrup) with ip and m being the 
azimuthal angle and quantum number, respectively, so 
two antennas placed on opposite sides of the billiard en- 
sure optimal coupling to all resonances. In the case of 
the square billiards, the antennas were placed offside any 
symmetry axes so that they couple to resonance states of 
all symmetry classes. It was checked that the results pre- 
sented here do not depend on the specific positions. The 
influence of waves reflected at the edges of the copper 
plates can be neglected. Details concerning the three dif- 
ferent dielectric billiards are given in the corresponding 
sections. Since the indices of refraction of the dielectric 
plates are only known with an uncertainty of a few per- 



cent, the precise values of n of the different plates were 
deduced from the length spectra by adjusting n such that 
the positions of the peaks in the length spectra match the 
lengths of the corresponding POs [53] , as will be shown 
below. The quality factors Qj = Re(/j)/rj of the mea- 
sured resonances are always smaller than those expected 
theoretically because there are Ohmic losses in the copper 
plates and the antennas and absorption in the dielectric 
material in addition to the pure radiation losses. 



IV. CIRCULAR TEFLON RESONATOR 

The first resonator investigated is a circular disk made 
of Teflon (Griinberg Kunststoffe GmbH). The resonator 
(called the Teflon circle in the following) has a radius of 
R — 274.9 mm and a thickness oi d = 5.0 mm. Its in- 
dex of refraction deduced from the length spectrum is 
n = 1.419 ± 0.001. Therefore, the critical angle for total 
internal reflection (TIR) is acrit — 44.8°. The dipole an- 
tennas were placed along the diameter of the disk on op- 
posite sides [see Fig. [l(b)[ . A measured frequency spec- 
trum is shown in Fig. [2j with a frequency of 10 GHz cor- 
responding to kR = 57.6. The spectrum features several 
families of almost equidistant, sharp resonances. These 
can be labeled with azimuthal and radial quantum num- 
bers (m, rir), and each family consists of resonances with 
the same rir and different m. The quality factors of the 
resonances are typically Q — 1000 — 5000. Since the 
radiation losses |Im(/m,„j,)| increase with increasing n^, 
only modes with small are observed in the measured 
spectrum. Modes with higher rir are visible at higher 
frequencies (compare, e.g., (18)). An angular momentum 
of hm can be attributed to each resonance. In the ray- 
picture, a trajectory with this angular momentum has an 
angle of incidence a with respect to the surface normal 
given by [23] 
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These trajectories are in general not POs. Two examples 
of trajectories associated with a resonance with Ur ~ I 
and with one with rir = 2 are shown as insets in Fig. [2] 
These show that a larger rir corresponds to a smaller 
angle of incidence a, because for a given resonance fre- 
quency [i.e., Re(fcj)], a mode with higher radial quantum 
number rir has a smaller azimuthal quantum number to, 
and a smaller angle of incidence results in larger radia- 
tion losses. The trajectories that are located close to 
the boundary of the circle have a large caustic (dashed 
inner circle) and are therefore called whispering gallery 
modes (WGMs). Indeed, all the long-lived resonances 
in the measured spectrum are of the WGM type. The 
bars above Fig. [2|indicate the frequency regimes in which 
modes with different radial quantum numbers can be ob- 
served. Most of the observed resonances have ri^ = 1 — 3 
and only a few have ri^ > 4. 
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FIG. 2. Frequency spectrum of the Teflon circle. The modulus of the transmission amplitude is plotted with respect to 
the frequency /. The classical trajectories associated with two resonances with Ur = 1 and rir — 2 are shown as insets, and a 
denotes the angle of incidence with respect to the surface normal of these trajectories. The bars above the graph denote the 
frequency regimes in which modes with certain radial quantum numbers Ur are observed. 
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FIG. 3. The integrated resonance density A'"(/) for the Teflon 
circle. The solid line is deduced from the measured frequency 
spectrum in Fig. [2j the dashed line is deduced from Weyl's 
law. All in all, 716 resonances were identifled in the spectrum 
up to 20.5 GHz, which is only about 10 % of the total number 
of resonances. 



In Fig. [3]the integrated resonance density N{f), which 
denotes the number of resonances below a given fre- 
quency /, is shown. The resonances observed in the 
measured spectrum (solid line) were counted twice be- 
cause all modes with m > are twofold degenerate. The 
dashed line is obtained from Weyl's law [Eq. Only 
modes up to 20.5 GHz < /2D are considered in the fol- 
lowing. The comparison of Weyl's law and the integrated 
resonance density in Fig. [3] illustrates that the long-lived 
resonances yield only a small part of the whole spectrum. 
The length spectrum for the Teflon circle is shown in 



Fig. |4j The full line is the length spectrum obtained 
from the measured frequency spectrum via Eq. ([s]), the 
dashed line corresponds to the complete spectrum calcu- 
lated by solving the Helmholtz equation for the dielectric 
circle , and the dotted line is the semiclassical expres- 
sion |psci(^)| from Eq. Q. All Fourier transforms were 
performed using the Welch-function as window function 
to smooth the resulting curves [5S] . The positions of the 
peaks of the experimental length spectrum depend sensi- 
tively on the index of refraction n used in the FT [cf. the 
term exp {—ikn£) in Eq. ([s])]. The positions of the peaks 
only coincide with the lengths of the POs if the correct 
value of n is used, and thus n is determined by matching 
these. The lengths of the different POs and the circum- 
ference are indicated by the arrows. The POs in the circle 
billiard have polygon and star shapes and are character- 
ized by their periods and rotation numbers (g,?]), where 
q is the number of reflections at the boundary and rj is 
the number of turns around the center, e.g., the (4, 1) 
orbit is a square and the (5, 2) orbit a pentagram. For 
the lengths considered in Fig. |4j only polygonal POs ex- 
ist (r? = 1). The orbits with q > 8 arc not indicated in 
the figure because their amplitudes Bpo decrease rapidly 
with increasing q (see Appendix |A]) . The semiclassical 
expression (dotted line) and the calculated length spec- 
trum (dashed line) agree very well except for the case 
of the square orbit. Its angle of incidence a = 45° is 
close to the critical angle for TIR. Thus, further cor- 
rections for Rpo must be taken into account |13j . No 
discernible peaks are visible for orbits with angles of in- 
cidence smaller than the critical angle. The overall shape 
of the experimental length spectrum is reproduced by the 
calculated length spectrum, but it has smaller peak am- 
plitudes, as to be expected due to the large number of 
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FIG. 4. Length spectrum for the Teflon circle. The full line 
results from the measured spectrum, the dashed line results 
from a complete, calculated spectrum, and the dotted line 
shows the semiclassical expression |psci(^)|. The arrows indi- 
cate the lengths of the depicted POs and of the circumfer- 
ence 27r_R of the circle. The semiclassical expression and the 
calculated length spectrum agree well except for the square 
orbit. The experimental length spectrum has smaller ampli- 
tudes than the calculated one. However, the deviations are 
smaller for the higher-order polygon orbits. 



missing resonances. In fact, in some cases the amplitudes 
of the peaks in the experimental length spectrum are as 
large as 80% of those in the calculated one. Thus, indeed 
the 10% most long-lived modes suffice to reproduce most 
of the peaks expected semiclassically. Interestingly, the 
agreement between the experimental and the calculated 
length spectra is better for the higher-order polygon or- 
bits such as the hexagon and the heptagon, and worst 
for the square orbit. This might be explained by the fact 
that the experimental spectrum consists only of WGMs 
(see Fig. |2|. 

To achieve a better understanding of the correspon- 
dence between different families of resonances and the 
different POs in the length spectrum, the measured spec- 
trum was divided into subspectra with radial quantum 
numbers rir = 1,2,3 and > 4. The radial quantum 
numbers were identified by comparison with the calcu- 
lated spectrum and by following the different series in 
the spectrum. First, only modes with = 1 are consid- 
ered, which form a family of almost equidistant resonan- 
ces. Note that spectra containing only one such family 
are often encountered in microlaser and microcavity ap- 
plications (e.g., [IS1I2S])- The length spectrum for the 
modes with = 1 is depicted in Fig. [5] Only reso- 
nances up to a certain frequency /max were considered 
for the three different curves: the solid line shows the 
length spectrum for all modes with = 1, the dashed 
line shows that for modes up to 15 GHz, and the dotted 
line shows that for modes up to 10 GHz. Each curve has 
only one peak, and, as expected, the position of the peak 
is close to the lengths of the high-order polygons, but the 
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FIG. 5. Length spectrum for a single resonance family. Only 
modes with radial quantum number = 1 were considered 
and only up to a certain frequency /max (solid line: 20.5 GHz; 
dashed line: 15 GHz; dotted line: 10 GHz). The arrows de- 
note the lengths of some POs and of the circumference 2t:R of 
the circle. Apparently, the position of the maximum of \p{i)\ 
only depends on the frequency /max, and is not related to any 
PO. 



position of the peak does not correspond to any partic- 
ular PO, instead it only depends on /max- This can be 
understood by analyzing the frequency spectrum in more 
detail. The resonance frequencies respectively wave num- 
bers for TM-modes with small Ur can be approximated 
as 
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where Xj is the modulus of the jth zero of the Airy func- 
tion Ai(a;) [27] . Therefore, the resonance spacing between 
modes with the same radial quantum number is 
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and the smallest resonance spacing (AA:)min is deter- 
mined by the highest azimuthal quantum number m or 
resonance frequency /max, respectively. Since the reso- 
nance spacing is almost constant, only one peak is ex- 
pected in the FT of the spectrum, and its position is 
related to (Afc)min- These observations are also made for 
the other two measured families [see Fig. 6(a)| . 

As a next step, we consider the length spectrum taking 
into account several families of resonances in Fig. |6(b)[ 
The combination of just two families (solid length spec- 
trum) is enough to obtain multiple peaks, whose posi- 
tions are indeed very close to the lengths of the POs. 
This effect can be explained as an interference between 
the different families (see Appendix [B|) . Moreover, the 
resonances with rir = 1 and 2 contribute mainly to the 
high-order polygons such as the hexagon and the hep- 
tagon, whose lengths are in the region where the corre- 
sponding single-family length spectra [full and dashed 
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FIG. 6. Length spectra for different resonance families, (a) 
Only one radial quantum number taken into account each 
(solid line: rir = 1; dashed line: rir = 2; dotted line: Ur — 3). 
(b) Several families combined. The solid line is the length 
spectrum for the families with radial quantum numbers = 
1 and 2 combined, the dashed line is for rir = 1,2, and 3 
combined, and the dotted line is for all resonances (identical 
with the solid line in Fig. |4]). The peaks of the length spectra 



(b) originate from the interference between the different 



resonance families, and the main contributions of the different 
families are focused on different orbits. 



lines in Fig. |6(a)| are maximaL The resonances with 
rir = 3 [added to obtain the dashed length spectrum 
in Fig. |6(b)| mainly contribute to the pentagon orbit, 
which again lies in the length region where the corre- 
sponding single-family length spectrum is maximal, and 
the few remaining resonances (added to obtain the dot- 
ted length spectrum) provide contributions to the pen- 
tagon and square orbits. So there is indeed a connection 
between the trajectories to which the different resonan- 
ces correspond and the POs in the length spectrum to 
which they contribute, and the deviations between the 
peak amplitudes of the measured and calculated length 
spectra are larger for the square and pentagon orbits (see 



Fig. |4]) because the most long-lived states correspond to 
the higher-order polygons. These results were also con- 
firmed by investigations with a complete, calculated spec- 
trum (not shown here). Additional losses due to, e.g., 
absorption in the Teflon material further reduce the am- 
plitudes of the experimental length spectrum, but this is 
only a secondary effect due to the generally high quality 
factors of the Teflon circle. 



V. SQUARE TEFLON RESONATOR 

The second resonator investigated is a square disk 
made of Teflon (called the Teflon square in the following) 
with side length a = 300.0 mm, thickness d = 5.1 mm 
and index of refraction n = 1.430 ± 0.001, which cor- 
responds to Qicrit = 44.4°. Its frequency spectrum is 
shown in Fig. [7| with a frequency of 10 GHz correspond- 
ing to ka = 62.9. The spectrum features a single fam- 
ily of broad, equidistant resonances atop an oscillating 
background. This background results from direct trans- 
mission between the two antennas. The quality factors of 
the resonances are in the range oi Q = 100 — 500, which 
is an order of magnitude lower than for the Teflon circle 
due to the larger radiation losses of the Teflon square. 
The resonances can be explained as superscarred states 
[TS] localized on the family of the diamond PO (inset of 
Fig. [7]). As the angle of incidence of the diamond orbit, 
Qpo — 45°, is very close to the critical angle, the radia- 
tion losses are large. The resonance frequencies can be 



approximated |15| up to O i 



as 



nL Re{k^) = mn + 2i ln[r(45°)] = ttto + 4(5 , (14) 

with L = \/2a being half the length of the diamond PO 
and TO being the longitudinal quantum number of the su- 
perscar. The phase S = arctan — Ijv?) is related to 
the Fresnel reflection coefficient via r(45°) = exp (— 2i5), 
and the term 4(5 accounts for the refiections at the dielec- 
tric boundaries of the square. The resonance frequen- 
cies computed from Eq. ( 14 ) are indicated by the arrows 
in Fig. [7] and show good agreement with the measured 
spectrum. Apparently, the observed resonances are su- 
perscarred states with longitudinal quantum numbers in 
the regime of to « 35 — 80 and first transverse excitation. 
States with higher transverse excitation are not observed. 
All in all, 49 resonances were counted up to 20.5 GHz, 
that is only 2% of the approximately 2220 resonances 
expected according to Weyl's law. 

In Fig. |8] the experimental length spectrum (solid line 
in the top graph) and the semiclassical expression |psci(^)| 
(bottom graph) for the Teflon square are shown. The 
arrows indicate the lengths of the different POs in the 
square, and dashed arrows are used for those POs not 
contained by TIR. The indices (jix^ riy) indicate half the 
number of bounces of the PO in x and y directions; see 
insets in Fig. [8) for examples. Details for the semiclas- 
sical amplitudes of the POs are given in Appendix [C] 
Only the (1, 1) orbit (diamond) and its repetitions are 
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FIG. 7. Frequency spectrum of the Teflon square. The series of equidistant resonances atop the slowly oscillating background 
are superscarred states localized on a family of POs including the diamond orbit (shown as inset). The arrows indicate the 
computed resonance frequencies of the superscarred states, and the corresponding quantum numbers m = 35 and m — 80 are 
indicated for two examples, respectively. 
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FIG. 8. Length spectrum of the Teflon square. The top graph shows the experimental length spectrum (solid line) and the 
length spectrum for a single family of superscars (dashed line), and the bottom graph the semiclassical expression. Note the 
different scales of top and bottom graphs. The arrows denote the lengths of the POs in the square (dashed arrows for orbits not 
contained by TIR), and the indices {nx,ny) denote the type of orbit. Two examples of POs are shown as insets [Fabry-Perot 
orbit (left) and diamond orbit (right)]. The experimental length spectrum features only the diamond orbit and its repetitions, 
which are the only POs confined by TIR. 



contained by TIR, and the experimental length spectrum 
only shows peaks at these lengths in agreement with the 
semiclassical prediction. Thus, the very simple structure 
of the frequency spectrum can be attributed to the single 
PO (plus its repetitions) contained by TIR. The trace for- 
mula also predicts a small peak at the length of the (1,0) 
orbit (Fabry-Perot orbit; left inset of Fig. |8]), which how- 
ever is not observed experimentally. The same applies for 
all other orbits not contained by TIR (indicated by the 
dashed arrows) . The peak amplitudes of the experimen- 
tal length spectrum are, however, less than 15% of the 
semiclassical prediction. Moreover, they decay exponen- 
tially with the number /i of repetitions, while the semi- 
classical amplitudes decay algebraically with 1/y/JI (see 
Appendix [C]). This is related to the facts that, first, only 



a single family of superscarred resonances is observed in 
the spectrum, i.e., the resonator essentially behaves like 
a one-dimensional system, and, second, that the resonan- 
ces have a finite lifetime. Starting with a single family of 
modes with resonance frequencies given by Eq. ( 14 1, one 
obtains 



2nL 



^ cos (2/^nLfc - 8/^(5) (15) 



as the semiclassical density of states 



|. The FT of 

/^fluc(^)' i-^-' length spectrum for a single family of 
superscars (ss) , is plotted as dashed line in the top graph 
in Fig. [8] and has a constant peak amplitude for all re- 
petitions /i of the diamond orbit. The agreement with 
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FIG. 9. (a) Ratio between the peak ampUtudes of the experimental and the superscar length spectra. The diamonds are 
the relative amplitudes for the indicated POs, and the full line is an exponential fit of Eq. ( |16[ ) whose slope corresponds to a 
resonance width of F = 56.1 MHz. (b) Experimental length spectrum (top graph) and superscar length spectrum multiplied 



with the exponential fit (bottom graph). The two curves show excellent agreement. 



the experimental length spectrum is better. However, 
Eq. ( 15 ) does not take into account that the system is 
open; in fact only the real part of k, i.e., the resonance 
frequencies, not the widths T are considered. The ra- 
tio of the peak amplitudes of the experimental length 
spectrum (full line) and the superscar length spectrum 
(dashed line) is shown in Fig. 9(a) with respect to the 
lengths ^po of the (repeated) diamond orbit. The solid 
line is a fit of the function 



Aq exp {—nipoTrV / c) . 



(16) 



with Aq = 1.09 and F = 56.1 MHz, which roughly 
matches the widths of the measured resonances. Indeed, 
Fig. |9(b)| shows that the experimental length spectrum 

(top graph) and the FT of p'l^^ik) multiplied with the 
fitted exponential decay corresponding to the finite life- 
times of the resonances (lower graph) agree with high 
precision. 

The example of the Teflon square demonstrates that a 
dielectric billiard with just a single dominant PO (the 
diamond orbit in this case) may act as an effectively one- 
dimensional system. Note that for a larger index of re- 
fraction or for higher frequencies, further families of su- 
perscarred states (i.e., with higher transverse excitation) 
might become visible in the frequency spectrum. Fur- 
thermore, the Teflon square is an example where a single 
family of resonances is directly related to one PO, which 
is in contrast to the case of the dielectric circle billiard 
(see Fig. [5] and corresponding text). The reason for this 
is that the resonances in the dielectric square are (su- 
per)scarred states, while those in the circle billiard are 
not. It appears that the relation of a single resonance 
family to POs depends on the specific case and needs 
careful analysis. 



VI. SQUARE ALUMINA RESONATOR 

The third resonator investigated is a square disk 
made of alumina (Morgan Advanced Ceramics) with side 
length a = 300.0 mm, index of refraction n = 3.050 ± 
0.008, and thickness d — 8.3 mm. The critical angle 
is acrit — 19.1°. A measured frequency spectrum of 
the resonator (called the alumina square in the follow- 
ing) is shown in Fig. [TO] with 5 GHz corresponding to 
ka — 31.4. The spectrum features many sharp resonan- 
ces, and 212 resonances were identified in the range of 
1.4 — 6.1 GHz (without taking into account possible de- 
generacies), compared to 1035 resonances expected ac- 
cording to Weyl's law. The quality factors are in the 
regime of Q = 200 — 2000. The radiation losses are gen- 
erally smaller than for the Teflon square due to the higher 
index of refraction, but the absorption losses in the alu- 
mina material are somewhat larger than in the Teflon 
material. In contrast to the case of the Teflon square, 
the spectrum shows no obvious structure of equidistant 
resonances, so the simple superscar model from Sec. |V] 
cannot be applied here because it takes into account only 
one family of periodic orbits with an angle of incidence 
of 45°. However, it can be generalized in a simple way 
|28j . We assume a ray with wave vector [k^^ ky) traveling 
in the square, where the x and y axes are parallel to the 
sides of the square, and the resonance condition after one 
round trip is 



r2(x) = l. 



«r2(7r/2-x) = 1, 



(17) 



with r(x) being the Fresnel reflection coefficient and x 
being the angle of incidence on the vertical sides. An 
approximate solution is 



k^a = m^TT + [r(x)] , 



kyd 



niy'K + i In [r(7r/2 — %)] 



(18) 
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FIG. 10. Frequency spectrum of the alumina square. In contrast to the cases of the Teflon circle and square, the structure of 
the spectrum seems irregular, and it cannot be divided into one or more families of equidistant resonances. The lower graph 
shows the part of the spectrum from 2 to 3 GHz. The arrows indicate the resonance frequencies computed according to a 
generalized superscar model, and the quantum numbers {mx,my) and the corresponding angles of incidence x are given for 
each. The computed resonance frequencies agree well with the measured spectrum in most cases. 



where {mx,my) are the x and y quantum numbers 
and the angle of incidence is approximated as x 

arctan (mj,/ma;). The wave number is k 



^kl + kl/n. 



This simple semiclassical approximation can be regarded 
as a generalized superscar model. A similar model was 
proposed in [291 . The case nix ~ Tny, X ~ 45° corre- 
sponds to the simple superscar model used in Sec. Mand 
resonance frequencies as in Eq. (14) are obtained. The 
resonance frequencies computed according to Eq. ( 18 ) 



are indicated by arrows in the magnified part of the spec- 
trum (lower graph in Fig. |10p . The quantum numbers 
are given together with x = arctan (ky/kx) ~ X, which is 
the angle of incidence on the vertical sides corresponding 
to the wave vector {kx,ky). In most cases, the measured 
and the computed resonance frequencies agree very well, 
and only in some cases they slightly deviate. Moreover, 
there are very few cases (e.g., at 2.7 GHz) where we can- 
not find a clear correspondence between a measured re- 
sonance and a computed mode. Only modes which are 
found in the measured spectrum are indicated. Modes 
with quantum numbers rux 7^ my are doubly degenerate. 
It should be noted that Eq. ( 18 ) predicts complex reso- 



nance frequencies [Im(fc) < 0] for x < Q^criti that is, losses 
due to refractive escape, and lossless modes [Im(fc) = 0] 
for X > c^crit. Note that in the latter case other ra- 
diative loss mechanisms, e.g., due to the corners [30] . 
are not accounted for. Accordingly, the model does not 
provide a complete description for the widths of these 
modes. All modes found in the measured spectrum have 
X > ckcrit- One would intuitively expect that the losses 
increase with x approaching ttcrit = 19.1° starting with 
X = 45° . Accordingly, the spectrum should be dominated 



by an equidistant series of resonances corresponding to 
X ~ 45°, as was observed in the case of the Teflon square. 
However, a closer inspection of the spectrum shows no 
clear correlation between x and the widths or amplitudes 
of the resonances. We attribute this effect to interac- 
tion between superscar and background states [31] . This 
would also explain why the computed resonance frequen- 
cies agree precisely with the measured ones only in some 
cases. Preliminary numerical computations confirm this 
interpretation of the effect, but it is not yet completely 
understood and will be the subject of further investiga- 
tions. In summary, the generalized superscar model ex- 
plains the general structure of the spectrum well, but not 
all of its details like, e.g., the resonance widths. 
The length spectrum for the alumina square is shown 
in Fig. |11[ The experimental length spectrum (upright 
graph in the top panel) features several peaks correspond- 
ing to POs, but most of the peaks hardly stand out above 
the noise level. The peak amplitudes of the experimental 
length spectrum are below 35% of those expected semi- 
classically (bottom panel), which is far less than for the 
Teflon circle, even though 20% of the total number of 
resonances was found in the case of the alumina square 
as compared to 10% in the case of the Teflon circle. A 
similar proportion between the number of resonances and 
the peak amplitudes was also found in numerical calcu- 
lations for the alumina square. The reason why a larger 
percentage of observed resonances in the case of the alu- 
mina square results in smaller peak amplitudes compared 
to the case of the Teflon circle is not understood, but we 
surmise that this is related to qualitative differences in 
the distribution of resonance widths for the circle and the 
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FIG. 11. Length spectrum of the alumina square. The upright graph in the top panel is the experimental length spectrum, 
the inverse graph is based on the generalized superscar model, and the bottom panel shows the semiclassical expression. Note 
the different scales of top and bottom panels. The lengths of the POs are indicated by arrows, with solid arrows for those POs 
confined by TIR and dashed arrows for those not confined. The indices {n^jUy) characterizing the confined POs are also given. 
There are several peaks corresponding to POs in the experimental length spectrum, but most of them hardly stand out above 
the noise level, and their peak amplitudes are much lower than semiclassically predicted. The inverse graph in the top panel 
is the length spectrum for a set of modes computed according to the generalized superscar model. Only modes with x ^ 28° 
were taken into account in this set, and no possible degeneracies. This length spectrum reproduces the experimental one well 
except for the cases of the (1, 1) and (2, 2) orbit. 



square billiard. No peaks corresponding to POs not con- 
fined by TIR (indicated by dashed arrows) can be iden- 
tified, as for the Teflon circle and square, respectively. 
For the observed ones the ratio of the experimental and 
the semiclassical amplitudes varies significantly. This is 
illustrated in Fig. [T2j where it is plotted with respect 
to the angle of incidence apo- Every trajectory in the 
square is characterized by two different angles of inci- 
dence (see Appendix [C|, of which the smaller one is used 
here. As in the case of the Teflon circle, even though they 
are still confined by TIR, POs with an angle of incidence 
close to the critical angle have significantly smaller rela- 
tive peak amplitudes compared to those with apo much 
larger than acrit- Finally, we compare the experimental 
length spectrum to the generalized superscar model in- 
troduced above. To this end, we calculate a spectrum 



according to Eq. (18), but omit all modes with angle x 
below a certain cutoff angle Xco, and include only one 
of the degenerate doublet modes since possible degenera- 
cies are not taken into account in the experimental length 
spectrum either. The cutoff angle is chosen as Xco = 28° 
to ensure that the number of resonances (215) obtained 



from 1.4 — 6.1 GHz is similar to that in the experiment 
(212). The lower part of Fig. [TO] shows that indeed al- 
most all measured resonances correspond to a x > Xco- 
Note also that the relative amplitudes in Fig. [12] start 
to decrease at a similar angle. Since the generalized su- 
perscar model does not provide a description of the re- 
sonance widths in the cases considered here, we assume 
a uniform resonance width of T = 8 MHz as obtained 
from an estimate of the typical resonance width found in 
the measured spectrum. The resulting length spectrum 
is shown as an inverse graph in the top panel of Fig. |11| 
It agrees well with the experimental length spectrum ex- 
cept for the cases of the (1, 1) and (2, 2) orbits, for which 
the amplitudes are significantly larger than those of the 
experimental length spectrum. These deviations cannot 
be explained yet. Nevertheless, even though the gener- 
alized superscar model cannot reproduce all details of 
the frequency and length spectra of the alumina square, 
most of the experimental findings are well described by 
this simple semiclassical approximation if taking into ac- 
count only the resonances which are observed in the ex- 
periment. It should be noted that the semiclassical trace 
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FIG. 12. The ratio between the peak amphtudes of the ex- 
perimental and semiclassical length spectra with respect to 
the angle of incidence of the corresponding PO. The indices 
(nxjUy) of the POs are indicated. The vertical line signifies 
the critical angle for TIR. The relative amplitudes for POs 
close to the critical angle are smaller compared to those of 
other POs. 

formula [Eq. Q] for the dielectric square can be deduced 
from Eq. ([I8| when all of the modes are considered [55] . 
The case of the alumina square demonstrates again that 
the long-lived, observable resonances contribute mainly 
to those POs with angle of incidence far away from the 
critical angle, indicating that corrections to the semiclas- 
sical trace formula must be taken into account for POs 
close to it. Furthermore, the larger index of refraction 
and thus greater number of POs contained by TIR are 
reflected in a more complicated structure of the frequency 
spectrum. In conclusion, not only the shape of a dielec- 
tric billiard, but also its index of refraction, i.e., the cri- 
tical angle for TIR, drastically influences the structure of 
its frequency and length spectra. 



VII. CONCLUSIONS 

We have measured the frequency spectra of three 2D 
dielectric microwave resonators with different shapes and 
indices of refraction. The length spectra were determined 
from the measured data and compared to the predictions 
of the trace formula proposed in [l3]. Even though the 
number of resonances identified in the measured spec- 
tra is only a fraction between 2% and 20% of the to- 
tal number of resonances, the experimental length spec- 
tra show peaks connected with the classical POs as pre- 
dicted, but the amplitudes of the peaks are smaller than 
the semiclassically expected ones. The differences be- 
tween experimental and semiclassical peak amplitudes 
depend strongly on the angles of incidence of the PO: 
POs which are not contained by TIR are not observed 
in the experimental length spectra, even if a finite am- 
plitude is predicted semiclassically. This was the case 
for all of the passive resonators investigated here, but 



there are indications that they may appear in the length 
spectra of strongly pumped microlasers [32j . Further- 
more, even POs contained by TIR have much smaller 
amplitudes if their angles of incidence are close to the 
critical angle. Apparently, there are two different effects. 
First, the strict semiclassical approximation derived in 
the limit /c — )■ oo, when applied for large but finite k, re- 
quires modifications which are especially important close 
to the critical angle. From the soluble case of the disk it 
follows that the reflection coefficient for such scattering is 
noticeably smaller than that predicted by the usual Fres- 
nel formulas [121 132- Second, it appears that long-lived 
resonances that can be observed experimentally mainly 
correspond to orbits which are not only confined by TIR, 
but also have an angle of incidence much larger than 
the critical angle. The examples of the Teflon and alu- 
mina square billiards further demonstrate that there is a 
connection between the number of POs contributing to 
the length spectrum and the complexity of the frequency 
spectrum. Therefore, it could prove useful for the de- 
sign of resonators to consider the type and number of 
POs confined for a certain billiard geometry and index of 
refraction. In general, the systematics of observed, long- 
lived and not observed, short-lived resonances must be 
taken into account for a full understanding of the exper- 
imental length spectra. Special care must be taken in 
the case of spectra with only one or two families of reso- 
nances. In general, there is no direct relation to certain 
POs. Only the case of (super)scarred states may be an 
exception to this rule, as the case of the Teflon square 
billiard demonstrates. 

In conclusion, we have demonstrated that the length 
spectra of dielectric resonators may be described by the 
trace formula, offering another tool to understand the 
correspondence between the ray- and wave-pictures of 
these devices. An advantage of the trace formula is that 
it only needs the spectrum, but not the field distributions 
as an input. On the other hand, the trace formula just 
gives information about a group of resonances in general, 
but not on the individual resonances, and at least a few 
dozens of resonances are needed in order to apply it. So 
far, the trace formula has merely been tested for passive 
systems with regular classical dynamics. Systems with 
chaotic or mixed dynamics or with an active medium 
(like microlasers) pose interesting future problems. 
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Appendix A: The periodic orbits in the dielectric 
circle billiard 



infinitely many orbits with length close to the circumfer- 
ence is vanishingly small. 



The POs in the circle billiard are characterized by their 
period and rotation numbers {q, rj), where q is the number 
of reflections at the boundary and rj is the number of 
turns around the center. Their lengths are 

fpo(g,r/) -2i?gsin(?77r/g), (Al) 

and the amplitudes and phases appearing in Eq. (|6| are 



Spo - ^\/2/po , With/po = I ^ ' 



1 , q^2rj^ 
otherwise, 



TT 

4 



q- +qarg [r{apo)] 



(A2) 



(A3) 



Here, r{a) is the Fresnel reflection coefficient for angle of 
incidence a with respect to the surface normal and elec- 
tric field perpendicular to the plane of incidence (TM 
polarization), and apo = Tr/2 — rjn/q is the angle of inci- 
dence of the PO. The area of the billiard covered by (the 
family of) the (9,77) orbit(s) is Apo — nR^ sin^ {rjn/q). 
These terms are identical to the case of a hard-walled 
circle billiard if r(apo) — —1 (see, e.g., [SUISj)- The am- 
plitude Bpo oc sin"^^^ (rjir/q), so Bpo — )■ for q/r] — !• 00. 
Therefore, the contribution to the length spectrum of the 



Appendix B: Length spectrum for two resonance 
families of the dielectric circle billiard 



Beginning with expression ( 12 ) for the resonance fre- 



quencies of the WGMs, a trace formula for the spectrum 
with just one or two resonance families can be derived. 
The (approximate) inverse of Eq. ( 12 ) is 



= Ak~ Bxnk^/^ + 



- 1 ' 



(Bl) 



with A = nR, B = {nR/2y^^, and Xj being the modulus 
of the jth zero of the Airy function. Following Chap. 3.2 
of we obtain 



Plr 



ik)=A(^A~^Bxr.^k-^/^'^ 



x^cos 2-K'q{Ak - Bxn.k^'^) (B2) 

»7=1 

as the fluctuating part of the density of states for a sin- 
gle resonance family. A factor of 2 was introduced to ac- 
count for the degeneracy of the resonances and the last 
term in Eq. (Bl ) was neglected. Adding up two of these 



terms for different radial quantum numbers and applying 
a trigonometric addition theorem results in 



1 



r;— 1 



2TT7^{Ak - Bxk^/'^)] cos [TrryBAxfc^/^ 



(B3) 



for two resonance families (e.g., = 1 and 2) with x = if its argument equals ^tt, i.e., for 
{x2 + xi)/2 and Aa; — X2 — xi. The first cosine term 
is a fast oscillating function, while the other terms are 
only slowly varying in fc, and the FT of p2n,. is computed 
using the stationary phase approximation. The result is 



X f TjAx 

6 



(B5) 



(7=1 



{2TrT]A - £n)^ 



TrriBAx\ 



' ^TTT/Bx 

2'KriA - in 



(B4) 



The details on the stationary phase (^^ are omitted. 
It should be noted that this expression (for each 77 
individually) is only valid if the stationary point lies 
within the integration interval, which is the case for 

£ < 2Trr]R[l~' 7^^{nR/2y^^xkmiL^]- For the same reason 
the positions of the maxima in Fig. [5] depend on /max- 
The modulus of the cosine term in Eq. (|B4| is maximal 



The numerical values of £max(97 v) close to the lengths 
£po{q, if) of ttic POs with (g, rf) for large (7/77, as shown in 
Table m which explains why the length spectrum for just 
two families of resonances in Fig. |6(b)| features peaks at 
the lengths of these POs. Indeed, it can be shown that 



ipoiq/v 00) = 27ri?77 



(B6) 



and that the zeros Xj of the Airy function fulfill 
x{Ax)'^ « TT^ for X = {xj+i + Xj)/2 and Ax — Xj+i — xj. 
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TABLE I. Comparison of the lengths l^o of the periodic or- 
bits in t he Teflon circle with the length £max according to 
Eq. (B5l. The values are computed for R = 274.9 mm and 
77 = 1. The agreement is very good for large g, and reaso- 
nable even for smaller q. 



q 


^po(g,r? = 1) (m) 


^maxlg,?? = 1) (m) 


A£ (m) 


3 


1.4284 


1.4126 


0.0158 


4 


1.5551 


1.5502 


0.0049 


5 


1.6158 


1.6140 


0.0018 


6 


1.6494 


1.6486 


0.0008 


7 


1.6698 


1.6694 


0.0004 


8 


1.6832 


1.6830 


0.0002 


9 


1.6924 


1.6923 


0.0001 


10 


1.6990 


1.6989 


0.0001 



Appendix C: The periodic orbits in the dielectric 
square billiard 

A PO traversing the square billiard 2nx i'2ny) times in 
the X (y) direction is denoted by the two indices (rix, riy). 
The length of the PO is 



where a denotes the side length of the square. Since all 
families of POs cover the whole billiard area, the ampli- 
tudes and phases entering Eq. ^ are 



Bpo — 



^ no 



,3/2 



ni + ni 



(C2) 



with Fpo = 2 in the case of diamond and Fabry-Perot 
orbits and Fpo = 4 for all other POs [35, and 



4 



arg (i?po) . 



(C3) 



ipoinx.riy) = 'ia^nl + , (CI) arctan (nj^/n^.) and a'^^ = tt/2 



Like for all regular billiards, Bpo oc and es- 

pecially Bpo oc 1/ where fj, is the number of 
repetitions of a PO. The total Fresnel coefficient is 
Rpo — r'^"^(apo)r-^"!'(Q!po), and the angles of inci- 
dence on the vertical and horizontal sidewalls are 

Qfpo, respectively. 
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